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EXECUTIVE  SUMMARY 


Fourier  transform  infrared  emission  spectroscopy  was  recently  implemented  for  in  operando 
studies  of  high  temperature  solid  oxide  fuel  cells  (SOFCs).  The  method  is  intended  to  detect  and 
spectrally  resolve  selective  emission  from  molecular  species  on  or  above  the  anode,  i.e.  for  hot  gases  over 
the  anode.  However,  the  substrate  and  other  broad-band  emission  sources  (e.g.  furnace  walls)  contribute 
to  the  signals  observed  and  complicate  their  analysis  and  interpretation.  These  effects  must  be  accounted 
for  in  order  to  estimate  the  sensitivity  to  selective  emission  and  to  segregate  the  molecular  components  of 
the  emission.  This  report  describes  the  effects  of  background  surface  emissivity  on  the  observed  spectral 
radiant  intensity  from  a  hot  gas  located  above  the  background  surface.  An  analytic  model  is  developed 
and  confirmed  numerically.  It  demonstrates  that  the  emitted  background  radiation  is  coupled  to  gas-phase 
molecular  emissions  in  a  nonlinear  way  and  cannot  be  completely  removed  from  the  molecular  gas 
spectrum  by  direct  subtraction  of  the  background  spectrum.  The  emission  spectrum  also  depends  upon  the 
experimental  geometry  through  the  radiative  form  factors. 
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The  Effect  of  Substrate  Emissivity  on  the  Spectral  Emission  from  a 

Hot-Gas  Overlayer 


INTRODUCTION 

Recently,  in  operando  optical  investigations  of  high  temperature  solid  oxide  electrochemical  cells 
(SOECs)  have  gained  significant  traction  with  several  research  groups  1-8  motivated  by  a  need  to  directly 
correlate  electrochemical  performances  of  SOECs  with  their  associated  and  complex  electrode  processes. 
These  techniques  (including  Raman  spectroscopy,  ambient  x-ray  photoelectron  spectroscopy,  near- 
infrared  thermal  imaging,  and  Fourier  transform  infrared  emission  spectroscopy)  have  demonstrated  a 
significant  advantage  over  traditional,  indirect  experimental  methods  by  providing  direct,  non-invasive 
information  in  real  time  along  with  high  temporal,  spatial,  and/or  molecular  resolution.  Fourier  transform 
infrared  emission  spectroscopy,  recently  implemented  for  in  operando  testing  of  solid  oxide  fuel  cell 
(SOFC)  anodes  by  Pomfret  and  coworkers,  4  has  become  a  novel  means  to  directly  measure  the  radiative 
emission  of  hot  gases  within  the  volume  over  an  anode  surface,  including  CFF,  CO2,  CO,  and  H2O  (see 
Figure  1).  In  these  measurements,  a  FTIR  spectrometer  was  adapted  by  removing  an  infrared  glow  source 
from  the  optical  path  and  aligning  in  its  place  the  anode  surface,  located  within  a  furnace  cavity.  Ideally, 
the  detector  will  only  see  the  infrared  emission  of  the  hot  anode  surface  and  any  gases  in  the  anode  head 
space,  though  in  reality  broad-band,  non-molecular  (and  frequency  dependent)  emission  sources  (e.g. 
furnace  walls)  contribute  to  the  majority  of  the  signal  reaching  the  detector.  These  must  be  subtracted 
from  each  spectrum  collected  during  an  experiment  to  leave  behind  the  difference  that  shows  only 
molecular  contributions. 

However,  recent  emission  experiments  have  indicated  that  for  the  FTIR  spectra  a  simple  background 
subtraction  may  not  account  for  all  background  contributions.  These  measurements  have  been  made  on  a 
variety  of  fuels  including  100%  CH4  and  mixtures  of  CEWAr  and  simulated  biogas  (50%  CH4  and  50% 
CO2,  balance  Ar)  at  various  operational  temperatures  which  represent  different  SOFC  operating 
conditions.  In  operando  IR  emission  spectra  collected  during  operation  with  both  fuels  described  above 
show  that  CO2  emission  is  significantly  attenuated  relative  to  other  gas  constituents  (Figure  2).  On  the 
other  hand,  ex  situ  IR  absorption  measurements  from  the  same  experiments  have  demonstrated  that  CO2 
absorption  is  optically  thick  and  very  intense  relative  to  other  gas  constituents  in  the  spectrum. 
Additionally,  attempts  to  calibrate  the  emission  intensity  of  CO2  with  concentration  in  a  high-temperature 
environment  indicated  a  nonlinear  response,  while  FT-IR  absorption  measurements  of  the  same  gases 
directed  into  a  flow  cell  showed  a  linear  response  as  expected.  Interestingly,  in  the  emission  experiments, 
as  the  concentration  of  CO2  increased,  the  normally  weak  CO2  combination  bands  at  3600  cm"1  and  3700 
cm"1  grew  in  intensity  while  the  asymmetric  stretch  at  2350  cm"1  remained  relatively  unchanged.  We 
expect  that  the  gas  composition  over  the  anode  headspace  is  comparable  to  the  exhaust  composition, 
despite  the  rather  large  apparent  concentration  differences  inferred  from  the  emission  and  absorption 
spectra.  To  better  understand  these  differences,  a  model  has  been  implemented  to  describe  how  the 
frequency-dependent  substrate  emission  could  affect  the  emission  of  a  hot  gas.  This  report  describes  a 
model  where  the  background  radiation  is  coupled  to  molecular  emission,  and  cannot  be  subtracted 
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linearly.  Furthermore,  the  geometry  of  the  test  apparatus  greatly  influences  any  background  radiation 
through  the  radiative  form  factor  and  reflections  from  the  anode  surface. 

A  similar  analysis  had  been  developed  previously  to  describe  the  conditions  necessary  for  selective 
emission  from  a  heated  solid  coating  on  a  substrate.  In  an  earlier  NRL  memorandum  report, 9  this 
problem  was  investigated  using  an  analytic  two-flux  radiation  model.  There  it  was  demonstrated  that  in 
order  to  have  selective  emission  from  the  coating,  the  emissivity  of  the  substrate  must  be  less  than  the 
emissivity  of  the  coating.  In  Figure  3,  taken  from  reference  9,  the  calculated  spectral  emission  of  a  1 10 
mil  thick  coating  of  Si02  on  steel  shows  selective  emission  since  Esioz  >  £steei  •  The  calculation 
assumed  that  substrate  and  the  SiC>2  coating  were  at  the  same  temperature  (1273  K). 
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Figure  1.  A  schematic  of  the  Fourier  Transform  Infrared  Emission  SOFC  test  assembly  is  shown  on  the  left  (a).  A 
more  detailed  view  of  the  headspace  between  the  anode  and  the  CaF2  window  can  be  seen  in  (b).  Modified  from 
original  figure  by  R.  J.  Kee. 


Effect  of  Substrate  Emissivity  on  the  Spectral  Emission  of  a  Hot-Gas  Overlayer 


3 


0.04  n 


O  0.03- 
c/) 

C/5 


E 

LU 


.  0.02  - 


C/5 

c 

<D 


In  Operando 


CO. 


BG,  no  load 


ft  »!V 


\j 


CO. 


2(ads) 


CH  polarized 


°-°H  ch4, 


a) 


no  load 

0  00  -W***^^ 

- 1 - i - . - < - r~  | — . — 

2500  2400 

2.0  n 


2300  220 

Wavenumber  (cm" 


2100 


o.o- 


BG,  no  load 


15'  CH.,  no  load 

0) 

O 

c 

05 

-8  1.0- 
O 
C/5 

< 

0.5- 


C02  hot  band  @  2270  cm 


1 — i 
2000 


Exhaust  Cell 
1 


_L.  u, 


JiJlyWV; 


Optically  thick 


17  torr  CO  in  Biogas:  0.60  torr  CO  in  CH  (28:1) 


3200 


3000 


2800 


2600 


2400 


2200 


2000 


b) 


Wavenumber  (cm'1) 


Figure  2.  FTIR  emission  spectra  collected  from  a  Ni/YSZ  anode-supported  SOFC  operating  with  simulated  biogas 
(BG),  CH4,  CO  and  CO2  at  800  °C  are  shown  in  (a)  and  corrected  with  a  linear  baseline  subtraction.  Optical  path 
length  was  approximately  5  mm.  During  the  same  trial,  anode  exhaust  was  fed  into  a  1 0  cm  gas  cell  and  the  absorbance 
was  measured  with  an  FTIR  spectrometer  (b). 


4  H.  D.  Ladouceur,  J.  D.  Kirtley,  D.  A.  Steinhurst ,  J.  C.  Owrutsky  and  S.  N.  Qadri 


Figure  3.  Selective  emission  from  SiC>2  coating  over  steel  substrate  at  1273  K.  Coating  is  1 10  mils  thick. 


RADIATIVE  TRANSFER  EQUATION 

The  background  contributions  observed  in  an  emission  experiment  could  arise  from  the  hot  walls 
of  a  furnace  via  reflection  or  a  hot  surface  oriented  normal  to  the  direction  of  observation  through  the 
emitting  gas,  e.g.  an  anode  in  an  active  solid  oxide  fuel  cell.  The  hot  emitting  gas  is  assumed  to  be  at  a 
spatially  uniform  given  temperature  T.  The  radiative  contributions  of  the  furnace  walls  and  other  surfaces 
within  the  observation  region  are  introduced  as  boundary  conditions  imposed  on  the  radiative  transfer 
equation10 

=  k(A)I(A,T)bb  -  k(A)I(A, T)  -  gscI(A,T)  +  ^  /  I(A,T)d>(A,  Si(  S)dfl, 

where  I  (A,  T)  denotes  the  spectral  radiant  intensity  (  radiative  energy  flow/time/area  normal  to  rays/solid 
angle/wavelength),  k(A)  denotes  the  spectral  absorption  coefficient,  I  (A,  T)BB  is  the  blackbody  Planck 
function,  osc  is  the  linear  scattering  coefficient  and  <f>(A,  §)  is  a  phase  function  associated  with  inward 

scattering  of  radiation  from  surrounding  control  volumes  into  the  direction  of  observation  about  the 
differential  solid  angle  dfL  This  term  increases  the  radiant  intensity  along  the  line  of  sight.  The  first  term 
on  the  right-hand  side  corresponds  to  spectral  emission  from  the  gas,  the  second  term  represents  spectral 
absorption  by  the  gas  and  the  third  term  accounts  for  scattering  losses  along  the  observation  path  s.  The 
scattering  processes  are  quite  small  in  the  absence  of  particulates  (Mie  scattering);  therefore  the  last  two 
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terms  will  be  neglected  in  this  analysis.  Assuming  the  hot  emitting  surface  is  located  at  s  =  0  ,  the  formal 
solution  to  the  radiative  transfer  equation  can  be  written  as 

I(A,  T,  s)  =  1(0,  A, Ts)exp(— k(A)s)  +  I(A,T)BB[1  -  exp(-K(A)s)]. 

The  factor  1(0,  A,  Ts)  denotes  the  spectral  intensity  from  the  substrate  located  at  s  =  0.  The 
transmittance  of  the  gas  x(A,  s)  over  a  distance  s  is  defined  from  the  Beer-Lambert  law 

x(A,s)  =  exp(-K(A)s). 

A0 

The  absorbance  A  is  defined  as  A  —  —  log(x(A,  s))  .  In  the  absence  of  scattering  or  reflections,  the  gas 
absorptivity  or  absorptance  a(A,  s)  is  given  by 

a(A,  s)  =  1  —  x(A,  s)  =  e(X,  s). 

By  Kirchhoff  s  law,  the  gas  emissivity  s(A,  s)  is  equated  to  the  absorptivity.  Using  these  relations,  we  see 
that  the  observed  spectral  line  intensity  I  (A,  s)  will  be  a  combination  of  both  the  transmitted  intensity  from 
the  surface  and  the  emitted  intensity  from  the  gas 

I  (A,  s)  =  I  (0,  A,  Ts)  exp  (— k(A)  s)  +  l(A,Tg)BB[s(A,s)], 

where  allowance  has  been  made  for  the  substrate  temperature  Ts  to  be  different  from  the  gas 
temperature  Tg. 

The  quantity  1(0,  X,  Ts)  represents  the  radiant  spectral  intensity  from  the  surface  at  a  given 
temperature  Ts  and  in  a  particular  direction.  The  calculated  spectral  intensity  depends  upon  what 
assumptions  are  made  about  the  substrate  surface  emissivity.  The  simplest  case  is  to  assume  that  the 
surface  is  a  blackbody  with  £  =  1  and  Ts  =  Tg.  Under  these  assumptions,  the  observed  spectral  intensity 
from  the  substrate  and  gas  overlayer  is  that  of  a  blackbody  at  the  common  temperature  T 

I(A, T, s)  =  I(A, T)BBexp(— k(A)s)  +  I(A,T)BB[1  -  exp(-K(A)s)]  =  I(A,T)BB  . 

In  the  second  case  let  the  surface  emissivity  be  that  of  a  grey  body  with  0  <  £s  (A)  <  1  and  assume  that 
there  is  no  incident  radiation  upon  the  surface  from  the  surroundings.  The  spectral  intensity  is  given  by 

I  (A,  s)  =  £s(A)l(A,Tg)BBexp(-K(A)s)  +  l(A,Tg)BB[e(A,s)]. 

This  equation  can  be  reduced  to  a  more  insightful  form 

I  (A,  s)  =  £s(A)l(A,Tg)BB(l  -  £(A,s))  +  l(A,Tg)BB[£(A,s)] 

=  (£sO)  +  e(^s)-£s  (A)  e  (A,  s))  I  (A,  Tg)  bb  . 
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This  result  shows  that  simple  linear  subtraction  of  the  substrate  background  from  the  spectrum  does  not 
completely  remove  the  background  due  to  the  cross-term  £s  (A)s(A,  s).  Upon  linear  subtraction  of  the 
background  contribution,  the  spectral  intensity  is  given  by 

lO,Tg)  =  e(X,Tg)l(A,Tg)BB[l  -  es(A)]. 

Note  that  the  “corrected”  spectral  intensity  is  given  as  the  product  of  two  factors.  The  first  factor 
e(\,  Tg)l(A,  Tg)BB  corresponds  to  the  molecular  emission  or  absorption  (Kirchhoff  s  law)  due  to  the  gas 

in  the  absence  of  a  substrate.  The  second  factor,  which  is  due  to  the  substrate,  reduces  the  spectral 
emission  intensity. 


MODEL  CALCULATIONS 

Figure  4  shows  an  idealized  emission  experiment.  A  nickel  substrate  is  placed  in  a  furnace  at 
1073  K  and  a  layer  of  10%  by  volume  CO2  is  located  above  the  substrate.  Using  the  HITRAN11  data  base 
and  a  commercial  code  for  calculating  molecular  spectra, 12  we  will  investigate  the  validity  of  the  above 
formulae. 
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Figure  4.  Geometry  for  idealized  emission  model. 


For  computational  purposes,  we  will  consider  the  CO2  antisymmetric  stretch  v3  centered  at  2349  cm'1 
(4.26  pm).  We  will  first  calculate  the  spectral  radiance  (W/nF/sr/cm'1)  emitted  by  the  gas  between  2200 
cm'1  and  2400  cm"1  in  the  absence  of  a  substrate  utilizing  the  HITRAN  line  parameters  and  the  Spectral 
Calculator  code.  The  results  of  this  calculation  are  shown  in  Figure  5. 
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Figure  5.  Calculated  CO2  spectral  radiance  without  substrate. 


The  calculated  in-band  radiance  (2200-2400  cm"1)  is  184.2  W/m2/sr.  Note  that  the  baseline  goes  to  zero  in 
the  wings.  Using  the  same  code,  the  blackbody  in-band  radiance  is  1389.2  W/m2/sr  as  shown  in  Figure  6. 
Therefore,  the  carbon  dioxide  emissivity  in  this  band  at  1073  K  is  given  by 


_  184.2 

£c°2  -  13g9 


0.1326 


2200 

Fri  Apr  10  1S.07“43  UTC  2015 


2250 


2300 

wavsnumter  (cm'1) 


2350 


P.  =  1013.25  mb 
T,  =  1073  K 
L  =  0.5  cm 


in-tand  radane* 
13S9.2  Wfmt sr 


2400 


Figure  6.  Calculated  in-band  Planck  radiance. 
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We  will  now  consider  the  effect  of  the  anode  surface  on  the  spectral  radiance.  The  anode  is  assumed  to  be 
made  of  nickel  with  an  emissivity13  £s  =  0.15.  The  calculated  spectral  radiance  is  shown  in  Figure  7. 


P  =  1013.25  mb 
T3'  =  1073  K 
L  j  =  0.5  cm 


in-band  radlanc* . 
364.979  Wftn  fer 


2200  2250  2300  2350  2400 

Fri  Apt  10  i  5.4Sj2$  UTC  20 1 5  waven  umbar  (cm ' 1 ) 

Figure  7.  Calculated  CO2  spectral  radiance  with  hot  nickel  substrate. 


The  calculated  in-band  radiance  is  365.0  W/m2/sr.  The  effective  emissivity  seff  of  the  anode-gas 
composite  system  based  upon  the  in-band  radiance  is  given  by 


£eff 


365 

1389.2 


0.263 


The  effective  emissivity  calculated  from  the  analytic  expression  (es  +  s(A,  s)-£s  e(X,  s))  equals  0.263 
which  demonstrates  the  validity  of  the  analytic  model  result.  Note  that  the  baseline  does  not  go  to  zero  in 
the  wings  due  to  the  substrate  emission. 


If  the  substrate  temperature  Ts  were  higher  than  the  gas  temperature  Tg,  the  observed  spectral 
radiance  would  show  an  absorption  feature  at  2349  cm"1.  If  the  transmittance  x(A,  s)  were  zero  in  the  2349 
cm"1  band  and  Ts^Tg,  the  observed  band  spectrum  would  be  that  of  a  blackbody  at  temperature  Tg.  If  the 
gas  temperature  exceeded  the  surface  temperature,  the  absorption  lines  would  be  inverted,  becoming 
emission  lines  above  the  cooler  blackbody  substrate. 
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Fraction  C02 


Figure  8.  Calculated  in-band  (2200-2400  cm"1)  radiance  as  a  function  of  CO2  volume  fraction.  The  black 
curve  corresponds  to  CO2  emission  without  a  substrate.  The  red  curve  includes  a  background  substrate  with 
an  emissivity  8=0.15.  The  gas  temperature  is  800°C  and  the  path  length  is  5  mm. 


Figure  8  shows  a  plot  of  the  calculated  in-band  (2200-2400  cm"1)  radiances  at  800°C  as  a  function  of  CO2 
volume  fraction  (f).  The  black  curve  was  calculated  using  the  Spectral  Calculator  code  over  a  5  mm  path 
length  with  a  total  pressure  of  1  atm.  The  red  curve  repeated  the  calculation  but  included  a  substrate  of 
emissivity  0.15  at  800°C.  Both  plots  show  a  linear  dependence  upon  concentration  for  small  volume 
fractions  of  CO2  (f  <  0.1).  As  the  CO2  volume  fraction  increases,  the  emission  spectra  show  a  distinct 
curvature.  This  behavior  is  consistent  with  the  Ladenburg-Reiche  function  for  the  curve  of  growth.10  At 
low  concentrations  the  behavior  is  linear.  At  higher  concentrations  the  asymptotic  behavior  follows  a 
square  root  function.  The  curve  of  growth  is  dependent  upon  the  line  strength,  the  concentration,  the  path 
length  and  the  temperature. 


EFFECTS  OF  SURROUNDING  SURFACES 

The  gas  emission  model  describe  above  is  a  one-dimensional  spatial  model  for  radiative  transport. 
The  radiation  sources  in  this  model  are  assumed  to  arise  solely  from  spectral  emission  of  the  substrate 
(anode)  and  the  overlayer  of  gas.  In  reality,  the  entire  system  would  be  embedded  in  a  furnace  with 
radiating  walls  and  windows.  The  spectral  emission  from  these  surfaces  can  contribute  to  the  observed 
radiance  by  spectral  reflection  from  the  substrate  surface.  The  radiation  incident  upon  the  anode  surface 


10  H.  D.  Ladouceur,  J.  D.  Kirtley,  D.  A.  Steinhurst,  J.  C.  Owrutsky  and  S.  N.  Qadri 


arrives  from  all  directions  and  is  termed  the  spectral  irradiation.  The  spectral  irradiation  is  denoted  by 
is  given  by 


r2n  rn/2 

=  I  I  I(A,  T,  0,  <^)cos0sin0d0dc|) 
J  o  Jo 


If  the  irradiation  comes  from  diffuse  surfaces,  I  (A,  T,  0,  cf>)  is  independent  of  0  and  cp  and  the  irradiation  is 
given  by 


G(A)i  =  tt  Ii(A,  T) 

where  the  subscript  denotes  the  i-th  enclosure  surface.  The  effects  of  reflection  can  be  included  in  the 
analysis  by  modifying  the  boundary  condition  for  1(0,  X,  Ts)  at  the  anode  surface  location  s  =  0. 

The  spectral  radiosity,  denoted  by  Js(0,  X,  Ts,  Tsurf)  accounts  for  both  the  emitted  spectral  energy  and  the 
reflected  spectral  energy  from  the  substrate  surface9 

Js(0,  A,Ts,Tsurf)  =  £sttIs(A, TS)BB  +  ps(A)^  G(A,Tsurf)i  F^s 

i 

where  ps(A)  is  the  substrate  spectral  reflection  coefficient,  G(A)i  is  the  incident  spectral  irradiation  on  the 
substrate  from  the  i-th  enclosure  surface  at  temperature  Tsurf  and  F^s  is  the  radiation  view  factor  defined 
below.  In  general  not  all  of  the  radiation  from  the  surrounding  surfaces  reaches  the  substrate  surface.  The 
fraction  of  emitted  radiation  from  the  surroundings  that  intercepts  the  substrate  surface  can  be  determined 
from  a  view  factor  or  configuration  factor  which  is  a  function  of  geometry.  The  view  factor  is 

defined  as  the  portion  of  radiation  leaving  surface  Si  and  intercepting  surface  S2,  that  is 

Radiation  leaving  Sx  and  intercepting  S2 
1->2  Total  Radiation  leaving  S± 

Catalogs  of  configuration  factors  are  available  for  diverse  geometries.14 

As  a  simple  application  of  a  configuration  factor,  consider  two  parallel  coaxial  circular  disks  of  radius  n 
and  r2  with  surface  areas  of  Ax  and  A2  separated  by  a  centerline  distance  h.  One  disk  (Fig.  9)  corresponds 
to  the  anode  surface  and  the  other  corresponds  to  the  calcium  fluoride  window  located  at  h  =  5  mm 
above  the  anode  surface.  The  configuration  factor  for  exchange  of  radiant  energy  from  to  A2  (anode  to 

window)  is  denoted  by  Fx^2  and  given  by  15 


The  configuration  factor  for  exchange  of  radiant  energy  from  A2  to  (window  to  anode)  is  denoted 
and  by  the  reciprocity  theorem16  is  given  by 
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Figure  9.  Anode -Window  Geometry. 


The  calculated  configuration  factors  for  this  geometry  along  with  the  spectral  properties  of  the  component 
materials  are  given  in  Table  1.  The  anode  is  treated  as  a  nickel  surface  with  an  emissivity  s  of  0.15  as 
described  above.  The  window  is  treated  as  a  grey  body  over  the  spectral  range  of  interest. 

Table  I  Anode -Window  Configuration  Factors  (h  =5  mm)  and  Surface  Properties 


Surface 

Radius  (cm) 

Area  (cm2) 

fm 

e(A) 

P(A) 

<A) 

Anode  (1) 

1.296 

5.277 

0.9082 

0.15 

0.85 

0 

Window  (2) 

2.002 

12.57 

0.3807 

0.03 

0.028 

0.94 

The  spectral  intensity  emitted  from  the  window  is  given  by  I  (A,  Twin)  =  £winI(A,  Twin)BB  where  £win  is 
the  spectral  emissivity.  The  spectral  radiosity  incident  on  the  anode  surface  is  given  by 

Gwin(^Twin)  —  ^win^K^  Twin)BBF2->i 

The  emissivity  of  a  3  mm  thick  calcium  fluoride  window  was  determined  from  an  absorption  spectrum. 
Figure  10  shows  the  derived  CaF2  window  emissivity.  The  feature  near  4.26  pm  (2349  cm"1)  is  due  to 
atmospheric  interference  from  carbon  dioxide  in  the  spectrometer.  Water  vapor  is  also  present  as  seen 
around  6.2  pm.  Using  the  values  from  Table  1,  and  the  expression 

£win^I(X  Twin)BBF2^lPs  (A) _  =  ^win^-^lPs  (A) 

Es*Is(A,Ts)bb  Ps rfwin)BBF2^i  £s  +  Ps (^-)^winp2^1 

for  the  fraction  of  incident  spectral  radiosity,  one  determines  that  the  window  contributes  6.1%  reflected 
radiation  to  the  observed  signal  from  the  anode. 
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Figure  10.  Emissivity  of  3  mm  thick  CaF2  window. 


The  calculation  of  radiation  from  the  cylindrical  side  wall  to  the  anode  surface  requires  the  numerical 
evaluation  of  the  configuration  factors  for  a  cylindrical  strip  (wall)  to  a  disk  (anode).  The  geometry  is 
shown  in  Figure  11.  The  outer  strip  is  5  mm  in  height  and  has  a  radius  of  2  cm.  Part  of  the  vertical  wall  is 
hidden  for  clarity.  The  configuration  factors  were  calculated  using  COMSOL’s  version  5.1  surface-to- 
surface  radiation  module  and  are  given  in  Table  2.  The  wall  is  AI2O3.  The  emissivity  of  the  AI2O3  was 
taken  from  reference  17. 


Figure  11.  Side  Wall  to  Anode  geometry 
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Table  2  Anode -Wall  Configuration  factors  and  Surface  Properties 


Surface 

Radius  (cm) 

Area  (cm2) 

g(A) 

P(A) 

x(A) 

Anode  (1) 

1.296 

5.277 

.09151 

.15 

.85 

0 

Wall  (2) 

2.002 

6.289 

.07681 

.40 

.60 

0 

Wall  (2— >2) 

.1172 

The  wall  to  anode  contribution  to  the  signal  is  14.8  %  for  the  assumed  geometry.  The  third  configuration 
factor  in  Table  2  arises  from  the  fact  that  the  outer  wall  “sees  itself’  in  the  radiative  transfer  process.  Note 
that  we  have  utilized  gray  body  properties  for  the  materials,  assumed  a  uniform  spatial  temperature  and 
treated  the  reflecting  surfaces  as  diffuse  reflectors.  As  the  wall’s  height  increases,  which  corresponds  to 
placing  the  anode  deeper  into  the  furnace  tube,  the  wall  to  anode  configuration  factor  will  increase  along 
with  the  background  radiation.  Given  the  spectral  properties  of  the  wall  and  anode,  the  spectral  emission 
intensity  could  be  calculated  over  a  range  of  wavelengths  using  the  same  algorithm. 

The  corrections  for  the  wall  effects  presented  here  assume  that  there  is  no  participating  (i.e.  absorbing  and 
emitting)  gas  located  within  the  enclosure.  The  corrections  would  be  applicable  to  surface  emissivity 
measurements  in  the  absence  of  a  participating  gas.  The  presence  of  a  gaseous  overlayer  modifies  the 
incident  radiosity  due  to  absorption  and  emission  along  the  propagation  path  and  complicates  the  analysis. 
There  are  several  methodologies  described  in  detail  by  Modest18  for  treating  this  type  of  problem.  These 
include  the  Monte  Carlo  method,  the  spherical  harmonic  method  and  the  discrete  ordinate  method.  This 
aspect  of  the  emission  will  be  the  subject  of  a  future  report. 


CONCLUSION 

An  analytic  model  for  the  spectral  emission  intensity  of  a  hot  gas  located  above  a  substrate  has 
been  developed.  The  substrate  contribution  to  the  observed  spectrum  is  nonlinear  since  the  effective 
spectral  emissivity  £eff  of  the  composite  system  is  given  by 

Eeff(A)  =  £S(A)  +  e(A,s)-£s  (A)e(A,  s) 

where  £S(A)  denotes  the  substrate  emissivity  and  s(A,  s)  is  the  spectral  emissivity  of  the  gas  layer  above 
the  substrate  surface.  The  product  £s  (A)£(A,  s)  is  on  the  order  of  one  percent  for  the  emissivities  used  in 
this  report.  Direct  subtraction  of  the  substrate  background  does  not  completely  remove  the  substrate 
interaction  with  the  gas  emission  spectrum.  Upon  linear  subtraction  of  the  background  contribution,  the 
“corrected”  spectral  intensity  is  given  by 

l(A,Tg)  =  £(X,Tg)[(X,Tg)BB[l-esW] 

The  first  two  factors  in  this  equation  correspond  to  the  selective  molecular  emission  of  the  gaseous 
overlayer,  which  is  linear  in  emitter  concentration.  The  third  factor  represents  the  effect  of  the  substrate 
emissivity  and  reduces  the  magnitude  of  the  “corrected”  spectral  intensity. 
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The  spectral  emissivity  of  a  gas  can  be  calculated  from  first  principles  utilizing  quantum 
mechanics  and  statistical  mechanics19.  In  this  report,  we  have  used  the  HITRAN11  data  base  in 
conjunction  with  the  Spectral  Calculator  software12  to  calculate  the  gas  emissivity  of  CO2  within  a 
spectral  band  centered  at  2349  cm'1  (4.26  pm). 

To  illustrate  the  validity  of  the  analytic  model,  the  in-band  radiance  for  the  carbon  dioxide 
antisymmetric  stretch  from  2200  cm'1  to  2400  cm"1  was  calculated  numerically  utilizing  the  commercial 
code  Spectral  Calculator  and  compared  with  the  analytic  model  prediction.  The  results  were  identical.  An 
overview  of  the  effects  of  the  furnace  walls  on  the  observed  emission  spectrum  has  been  presented.  The 
analysis  is  very  dependent  upon  the  detailed  experimental  configuration  and  material  properties  such  as 
emissivity  and  the  nature  of  the  emitting  surface  (diffuse  or  specular  reflector).  The  observed  spectral 
emission  can  depend  upon  the  location  of  the  system  within  the  furnace  since  emitted  wall  radiation  can 
interact  with  the  gas  and  the  substrate  through  the  radiative  configuration  factor  and  can  contribute  to  the 
observed  spectral  radiant  intensity.  The  reflective  contribution  is  on  the  order  of  15  %. 

The  contributions  from  the  surrounding  surfaces  to  the  spectral  intensity  described  here  do  not 
include  the  effects  of  a  participating  gas  layer  interacting  with  the  radiation  from  the  configuration  factor 
contributions.  Such  an  analysis  requires  the  numerical  techniques  described  by  Modest18  and  can  be 
implemented  in  COMSOL.  However,  the  analysis  is  applicable  to  corrections  for  surface  emissivity 
measurements. 
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